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1. Introduction 

The run of the Large Hadronic Collider at CERN with energies, at least an order of 
magnitude higher than the previous colliders, poses new challenges to the high energy 
physics community. One of them, faced by the experimentalists, is the handling of very 
high multiplicity events. The time needed to analyze such events grows more than linearly, 
when fitting procedures are employed, for instance - jet identification by clustering. 

On the other hand, as it will be shown, the higher the energies of the jets are, the 
measured accuracy of the momentum tensors' global variables increases, making them 
more suitable for analysis. 

Two quark-antiquark (qq) jets were first seen in e'^e~ events as two clusters of 
particles moving back to back in opposite directions, with the same absolute value of 
momentum. 

Three jets, with a gluon jet radiated from one of the quark jets, were observed later 
on, also in e~^e~ events [§, ^, ^. This is illustrated in Fig. ||. Similarly to two jets, the 
three jets appear in their centre of mass system with vanishing total momentum. 

For the study of jet events, a symmetric tensor of rank two, formed by the momenta of 
the particles, was suggested by |^. It was studied and modified later on by and 0, in or- 
der to obtain a tensor, free of mass singularities, and therefore computable by perturbative 
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QCD: 
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The summation over n is up-to the total number of particles N, p„ is the momentum of the 
particle n and the indices i and j indicate components of the momenta. The denominator 
serves for normahzation. 

The reference actually gives a gen- 
eral definition for such a symmetric mo- 
mentum tensor of any rank r: 



Qijk.. 



N 

E 

n=l 



Pn,iPn,jPn,k'''Pn,r 



N 

E \Pn 

n=l 



[1.2) 









- 



Figure 1: Schematic view of a 3-jet event 



Note that the normalization factor is iden- 
tical to that of rank 2, and the terms of the 
tensor are linear with respect to momenta. 

In addition to the important theoretical features, it can be easily shown that for any 
rank, the tensor remains invariant if one particle is considered as two collinear ones, or if 
two collinear particles appear as - one. Such errors occur frequently, when calorimeter data 
is used for neutral particles. 

Using the common convention that a repeated index means summation, one obtains 
that the trace of the tensor of rank two ( |1.1| ) is: 



Qii — Ai + A2 + A3 — 1 
where the A's are the eigenvalues. By writing them in decreasing order: 

Ai > A2 > A3 



(1.3) 



(1.4) 



and by defining an idealized jet with collinear particles only, meaning that a jet is expressed 
by a single momentum, one obtains that an event with two such jets has 



Ai = 1 A2 = A3 = 
and any idealized three jet event is planar: 

A3 = Ai > A2 > Ai + A2 



1 



(1.5) 



(1.6) 



Note though, that not all planar events correspond to three jets. 

The remnants of the beam particles in hadron-hadron collisions carry part of the 
energy and momentum, before disappearing inside the collision pipe. As a result, the two- 
jet events are not back to back, and the three jet events are not planar anymore. In order 
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to use the same tensors in the analysis, it is imperative to make the appropriate Lorentz 
boost before-hand, by requiring that the total momentum of the particles, originating from 
the event 's vertex, vanishes as a result of the transformation. 

For a more detailed study of three and more jet events, clustering algorithms were 
developed. A review article of algorithms for e^e~ events can be seen in Q. Algorithms 



for hadron-hadron events have been more recently developed, e.g |10|. A fitting procedure 
is used, which could be time consuming for high multiplicity events, as those obtained at 
the LHC. The present study could facilitate this procedure in the case of 3-jet events, by 
doing a priory selection. In addition it could yield some independent measurements for 
comparison. 

This paper is organized as follows. Section |2| deals with idealized jets and their prop- 
erties, analyzed by the momentum tensors of ranks two and three. Section ^ deals with 
more realistic jets, generated by Monte-Carlo, and concludes with suggested implementa- 
tions for real data. Section contains a summary and conclusions. 



2. Idealized jets 

An idealized jet, as defined in the introduction, has only collinear particles and is expressed 
by a single momentum. 

2.1 Tensor of rank two 



In the case of a planar event, the eigenvalues of the rank-2 tensor of momenta (1.1) fulfill 
the relations (|1.6|). Therefore, only one free parameter defines the global event shape, 
which is chosen here as: 

fi = XiX2 0</x<^ (2.1) 
Details of the physical meaning of the tensor's invariants can be found in 







S 

N 

\ 

N 

\ 

S 

s 

>-V-? 


-1 






























^ « Pi fi^ 







Figure 2: Idealized 3-jet event. 
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The simulation of an idealized 3-jet event is done here in the (x,y) plane by the momenta 
components of the jets: 



|pi(-pi,0)| > |p2 (P2Cosa,p2sina)| > |p3 (ps cos /3, -ps sin/3)| >0 (2.2) 

where the size of momenta are ordered for convenience. From the triangle of the momenta 
in Fig. m one obtains: 

\Pl\ \P2\ \P3\ „^ 

= = -■ = PO 2.3 

sm (a + p j smp sma 

where Pq is a free parameter of no importance and the angles fulfill: 

0<a</3<7 = 7r-a-/3 (2.4) 

The none-vanishing terms of the momentum tensor ( |1.1| ) corresponding to ( |2.2D are: 
DQii = D — sin a sin /3 (sin a + sin /3) 



DQ22 = sin a sin j3 (sin a + sin /3) 
DQ12 = sin a sin (3 (cos a — cos /3) 
D = sin a + sin /3 + sin (a + /?) 

After some trigonometric manipulations one obtains: 

sin a sin /? sin (q + /3) 



(2.5) 



// = A1A2 = 
with eigenvalues of the tensor: 



Qu Q12 
Q21 Q22 



sin a + sin /3 + sin (a + /?) 



(2.6) 



where from ( |2.6D can be rewritten also as: 

^ = cos|(cos^-cos0 where {^^^^"} (2-8) 

The domain of the angles a and /? is displayed in Fig. ^ together with some of the curves 
corresponding to constant values. The dots represent calculated values connected by 
straight lines. The following relations correspond to the boundaries: 

a = f3 cos/3 = = ^1 ^ |P2| = |P3| 

/3 = 7 cos/3 = ^ ^ = IP2I 



(2. 



The value of //, does not allow an unambiguous determination of the angles, except in the 
case of // = 0.25. The value of // = corresponds to a two-jet event, and therefore very 
small fj, values are not a good indication for a three-jet event. For a given ^ value, the 
angles a and f3 are constrained by the relation (|2.8D, and bounded to the values of (|2.9D. 



The /i values of planar events are bounded by (2.1), exactly as in the case of three 
ideal jets. In the latter case the maximal fi value corresponds to rotationally symmetric 
distribution of the jets. Let's define an idealized, rotationally symmetric planar event of 
N jets, as having identical absolute values of momenta, and each jet being separated by 
an angle of 27r/N from its neighbours. 
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Figure 3: Level lines of /i. See (|2.lD, ( [2.2| ) and ( p.4| ) for details 



An example of = 5 is illustrated in Fig. ^, 
and the momenta of the jets for any A > 3 are 
expressed by 



n 



[pocos(^) ,posin(2fi)] 



1, ...,A 



A > 3 



(2.10) 



In such a case, because of the symmetry, the 
eigenvalues are equal (Ai = A2), and therefore 
from (2.1) and ([L.6|), /x becomes maximal: 



A* 



(2.11) 




2k/5=72' 



X 



2n/5=72' 



This result can be obtained of course also by 
direct calculation, which is shown in App. in 
order to be compared later with tensor of rank Figure 4: Symmetric planar event of five jets 
three. 

2.2 Tensor of rank three 

From the definition of linear tensors ( |1.2| ), the tensor of rank three is: 



Qijk 



N 
n=l 

N 

E \Pn\ 
n=l 



(2.12) 



and is used in Q as an example of QCD calculation for 3-jet events. As far as we know, it 
has not yet been used as a tool for experimental analysis. 
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Contrary to the tensor of rank two, and due to the products of odd number of momenta, 
this tensor vanishes for any ideahzed two-jet event. This property, together with the tensor 
of rank two, can be used to enhance the signature of a two-jet event. 

From ( p. 12 ) and using the summation convention one obtains a vector: 



Vj — Qji 



N 

E 

n=l 



_ n=l 

~ N 

E \Pn\ 
n=l 



N 
n=l 



(2.13) 



that vanishes, since it is the total momentum of the particles in the centre of mass system. 
This yields three linear relations between the terms of the tensor: 



Qui + Q122 + <5l33 = 
Q2II + Q222 + <5233 = 
<5311 + Q322 + <5333 = 

and of course the absolute value of the vector also vanishes: 



(2.14) 



V,V, = 



One can also obtain a tensor of rank-2: 



(2.15) 



Rij — QimnQ 



]mn 



which yields: 



Rii 
R22 
R33 
R12 
R23 
R31 



Qiu + Q122 + Q133 + 
Q222 ~^ Q211 "I" Q233 ~^ 
Q333 + Q311 + Q322 + 
Q133Q233 — Q133Q211 
Q211Q311 — Q211Q322 

<3322<5l22 — Q322Q133 



2 (Q2II + Q123 + Q3I1) 
2 (Q122 + Q123 + Q322) 
2 (Q133 + ^123 + Q233) 

- <5l22Q233 — 2(5i23(5333 

- <5233Q311 — 2(5i23(5lll 

- <531lQl22 — 2(5l23Q222 



(2.16) 



(2.17) 



For planar events, assuming e.g. the (x, y) plane, any term of the Qijk tensor vanishes, 
if at least one of its indices equals 3. From ( p.l7| ) one obtains in such a case that there are 
only two non- vanishing terms of Rij: 



Rn — Qin + Q122 + 2Q211 



2 (Q?22 + Qlu) = 2 (Q?ii + Q 



222 J 



R 



'22 



(2.18) 



where the relations of ( p. 14 ) were used. If the eigenvalue of the tensor Rij are denoted by 
the letter ly and are ordered by decreasing values, in case of planar event one obtains: 

1^1 = 1^2 > 1^3 = (2.19) 
For a simulated idealized 3-jet event in the {x,y) plane, one obtains from ( |2.18D and ( 2.ig| ): 

Qlii + Q222 = (2.20) 
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where fi has exactly the same value as in the case of (i2.q). Using the eigenvalues one can 
summarize that for an idealized 3- jet event: 

u, = u, = 2X1X1 = 2,- I 
This is a genuine signature of a three-jet event, which was obtained by two entirely 



different mathematical mechanisms, while in general ( 2.20 ) does not hold. As an example, 
any planar rotationally symmetric event of > 3 jets fulfills A1A2 = /x = 0.25 (App. |A|) , 
but only for A^ = 3, ( |2.21 ) takes place, while for A^ > 3: vi = U2 = ^ (App. P). 



3. Real jets 

In a real event, due to the hadronization process, the momenta of the particles of a jet 
are close to the jet axis direction, but are not necessarily collinear with it. Consequently 
the planar property of a 3-jet event is broken, i.e. the smallest eigenvalues of momentum 
tensors, A3 and 1^3 ( |2.21 ), do not vanish any- more. 



In order to study this behaviour, 3-jet Monte-Carlo events were generated and particles 
produced by the hadronization process. The momenta and energies of the primary particles 
obtained from the hadronization were recorded in their common centre of mass system. 
They were used for the calculation of the momentum tensors, and for comparison with the 
idealized case. By definition, this system is equivalent to the common centre of mass of the 
jets. 



We have used version 6.408 [12| of the PYTHIA program to generate high-energy 
physics events. The given total energy, in the centre of mass system of the jets, was 
randomly distributed among the three jets. The most energetic jet was assigned to a quark 
jet, while the other two - to a quark and to a gluon jet randomly. The quark jets originated 
from n or d quarks. Events with a jet energy less than 2 GeV where not accepted, which 
resulted in the suppression of events with /i < 0.05, less pronounced at higher cm. energies. 

Some of the parameters, corresponding to the idealized events, can be obtained with 
a good approximation directly from the real events. Other require fitting and parameter- 
ization depending on the choice of the event generation mechanism, which is very simple 
here. For any other choice, the same procedure can be easily repeated, and the appropriate 
parameters obtained. 

3.1 Tensor of rank tw^o 

It is natural to start by asking how to obtain a measured value of ^ui, close enough to the 
idealized value of /x (^H]). Let's consider two definitions of /xi: /x'^ and fj,'( ( |3.1[ ). Since 
A3 ^ for real events, fi'i is smaller than the idealized /x, due to the trace invariance ([1.3|). 

/x; = AiA2 < M < A1A2 + A2A3 + A3A1 = ^ = /x'/ (3.1) 

The /x" is actually a tensor invariant, well studied in Q as the C parameter. For the real 
case, this value is bigger than that of the idealized case, due to the hadronization. For the 
idealized case, it also corresponds to /x (^). 
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The distributions of the mean A// = 
— fii values, where n corresponds to the 
ideahzed 3-jet events, and yUi to fi'i or to /i" 
( |3.1| ) , obtained from generated 40000 events 
for each one of the five cm. energy E, are 
shown in Fig. ^. Two different signs mark 
the data points. For ^u" (bottom), the values 
of A// are called hadronization corrections. 
The hadronization corrections of the event 
shape parameters, e.g. C (|3.1| ), have been 
studied first by [^] and [0] and developed 




40 80 120 160GeV 
_l I I I 



further by more authors. Reference |15] con- 
tains a review of these studies . The bottom 
distribution of Fig. ^ should behave approx- 
imately as 

Afi = fi- (A1A2 + A2A3 + A3A1) = --^ 

(3.2) 

where a is a positive constant. The two solid ^';^^p:^^ ^.^S^^ ^"^^^^^t^ ^^^'"^g^ 1^°*^ 
... ,11 1 ■ r ,1 distributions. The curves are explained in the 

curved lines express the boundaries of the 

^ text, 
data points, where the extreme values of a 

are used: 

0.75 < a < 0.92 GeV 



Figure 5: Distribution of the mean values of 
A/z (ordinate), as defined in the figure, vs the 
cm. energy of three jets (abscissa). The square- 



The top distribution of Fig. |5| behaves similarly, but with inverse sign, yielding: 

Afi = fi- A1A2 = -| 
0.74 <b< 0.82 GeV 



(3.3) 



(3.4) 

The values of a (|3.3D and b (|3.4[) are close to each other, and the graphical displays behave 
almost as mirror images. 

This behaviour can be explained by simple considerations. An idealized (i) 3-jet event 
in the {x,y) plane fulfills (^). 

A3(0 = Ai(i) > A2(i) >0 Ai(i) + A2(i) = 1 (3.5) 

After hadronization, vector components outside the plane appear, and one obtains: 

A3 = A3>0 ^ Ai = Ai(i)-Ai and A2 = A2(i) - A2 = 1 - Ai(i) - A2 (3.6) 

From the trace invariance ( |1.3D , one obtains: 

A3 = Ai + A2 (3.7) 

For given values of fi and E, the values of Ai and A2 could in principle differ one from the 
other, but by symmetry considerations we assume that their means are equal: 



Ai = A2 = A and A3 



2A 



(3.8) 
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(3.10) 



Assuming that A is much smaller than one, we can write that 

< A (3. 

yielding the first approximations of the mean values 

M = A1A2 = [Ai (i) - A] [A2 {i)-A] = fi-A [Ai (i) + A2 (i)] + A2 = ^ - A 
/i'/ = fi[ + A3 (Ai + A2) = /i'l + A3 (1 - A3) = /u - A + 2A (1 - 2A) = /i + A 

The A from ( p.lCl| ) is the mean value of the hadronization correction for 3-jet events, which 
can be measured experimentally by 

^ ^ ^ A3 (A1 + A2) ^ A3(i-A3) 

2 2 2 ^ ' ' 

An approximation of the idealized ^ value can be obtained experimentally by the mean 
value of /ii: 

+ 2AiA2 + A3(Ai + A2) 2AiA2 + A3(l-A3) 

Ml = 7, = 7, = 7^ I3.i2j 
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Figure 6: Distribution of the mean values of 
the hadronization corrections A (ordinate), 
defined in the figure, vs the cm. energy of 
three jets (abscissa). The curves represent 
the boundaries of the points. 
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Figure 7: Distribution of the mean value 
of A/^ (ordinate), defined in the figure, vs 
the cm. energy (abscissa) of three jets. The 
straight lines connecting the data points are 
just for guiding the eye. 



The hadronization corrections A (3.11), averaged over different energies , are plotted 
m Fig. 1. The curves that bound the plotted points are 



0.78 



< A = - < ^ (3.13) 

E - E - E ^ ' 

in agreement with ( |3.2D - (^^). Here they represent the values obtained from the event 
generator used. In the case of experimental data of 3-jet events, they should indicate the 
measured value of the hadronization corrections for fi 
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The arithmetic average of the upper and 
bottom distributions of Fig. ^, displayed by 
squared marks, corresponds to the A// = 
/X — /ii ( p. 12 ) distributions. The enlarged- 
scale distribution of the latter is given in 
Fig. |. 

From the same sample of events, the 
standard deviations of the distributions were 
calculated and plotted in Fig. ^. They ex- 
hibit values, more than one order of magni- 
tude larger than the corresponding A// val- 
ues of Fig. 0, therefore the choice of ( p. 12 ) is 
acceptable as an approximation to the ide- 
alized value. The solid curves from Fig. ^ 
bounding the distribution of the standard 
deviation are 



0.40 <d< 0.51 GeV 



(3.14) 
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Figure 8: Distribution of the standard devia- 
tion a of A/x (ordinate) defined in the figure vs 
the cm. energy (abscissa). The straight lines 
connecting the data points are just for guiding 
the eye. 



Fig. P and Fig. show scatter plots 
of the values of /U (ordinate) vs (|3.12|) 
(abscissa) obtained from tree-jet generated 
events, for two cm. energies: 80 and 200 

GeV. Note the depletion of events with at // < 0.05, mentioned [previously , and the nar 
rowing of the distribution with energy and with the increase of the ^ values . 
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Figure 9: /i (ordinate) vs /ii (abscissa) of Figure 10: /i (ordinate) vs /ii (abscissa) of 
3-jet events of 80 GeV cm. energy. 3-jet events of 200 GeV cm. energy. 
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Figure 11: Distribution of A/ii = /i — /ii ( 3.15 ) for 40GeV cm. energy 
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Figure 12: Standard deviation tr (ordinate) vs /ii ( |3.12 ) (abscissa) of 3-jet events of different cm. 
energies. The curves are fitted ellipses for the parameterization of a ( 3.161 ). 



Fig. 11 displays the distribution of 



(3.15) 
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obtained from 40000 generated events at 40 GeV cm. energy. Its mean value (A/i) = 
0.000771 represents the systematic deviation of ni ftom /U. The corresponding standard 
deviation a = 0.01012 is due to the hadronization process of the jets. 

In order to be able to use the standard 
deviation, one needs to express it by the mea- 
sured quantities: ( |3.12| ) and the cm. en- 
ergy E. For this purpose, the events of each 
one of the investigated energies (40, 80, 120, 
160, 200 GeV) were divided into four regions 
of m: 0.05 to 0.10; 0.10 to 0.15; 0.15 to 0.20; 
0.20 to 0.25. The region below 0.05 was not 
used, because of the depletion of events there. 
For each energy and region of /xi, the appro- 
priated standard deviation was calculated. 

The standard deviation a as function of 
the energy and /ii, obtained that way, is dis- 



played in Fig. 12. The solid lines represent 
elhpses, fitted by "Minuit" 0, yielding = 
3.1 for 15 degrees of freedom, and are ex- 
pressed by: 

rc = .25-/ii 
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Figure 13: The standard deviation a (ordi- 
nate) obtained by the statitical distribution, 
and by the parametrized expression, vs the 
cm. energy (abscissa) 



1 = ao[E{E + ai) + 02] 
(3.16) 

where the cm. energy of the event E is in GeV, and the numerical values of the parameters 
are in the appropriate units: c = 1.812 ± .347; e = -25.74 ± 2.85; oq = -.00001485 ± 
.00000486; ai = -355.3 ± 27.4; 02 = 4864 ± 2839 . Ehipses were chosen just for practical 
reasons, in order to obtain a simple formula for interpolating the fii and E values. The fixed 
point of fj = at ^1 = 0.25 for all energies is backed up by the trend of the data points, 
consistent with being a stationary point corresponding to the maximal ni. The obtained 



a values from ( 3.16 ) are due to the hadronization process alone, and the experimentalist 
should take care of adding the errors of measurement. E = 40 GeV is used as the lower 
limit for the present analysis and any extrapolation below this energy could be erroneous. 
This limit was chosen, since for lower energies, the standard deviation increases sharply, 
and the use of fii becomes non-practical. On the other hand an extrapolation to higher 
than 200 GeV energies is expected to be safe, fii < 0.05 was not used for the fits, and it 
is advisable to make at least such cut in the data, in order to avoid two-jet events. 

A test of the standard deviation function (|3.16 ) was done by calculating o" for A'' = 
40000 generated events at each given energy, by use of 



N 
n=l 

N 



(3.17) 
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where (T„ corresponds to the value of a ( 3.16 ) for the n*'^ event, rather than from the 
deviation of the statistical distribution, presented previously at Fig. |^. The results, super- 
imposed on the statistical calculation, are shown at Fig. |l^ for comparison. Both exhibit 
a pattern close to each other, with deviations of less than 10%. This is satisfactory in view 
of the two different procedures, giving us confidence in the parameterization (3.16). 



3.2 Tensor of rank three 

It was shown that the use of the rank-3 tensor, for idealized jets, supplies together with 
the tensor of rank-2, a signature of a 3-jet event. In such a case the eigenvalues Ui of the 
matrix Rij (2. IE) fulfill ( |2.21| ), which is reproduced here for convenience: 



1^1 

^3 



V2 

As 



2Af A2 

: 



2/i^ 



(3.18) 



where Aj are the eigenvalues of the tensor of rank-2 Qij (1.1). Real events do not fulfill 
vi = V2 anymore, but they remain close to each other, and vj, does not vanish. 




Figure 14: /i (ordinate) vs /i' (abscissa) of Figure 15: /i (ordinate) vs p! (abscissa) of 
3-jet events with 80 GeV cm. energy. 3-jet events with 200 GeV cm. energy. 

As a first approximation of the value of /U obtained from the tensor of rank-3, we take 
according to ( |3.18| ): 



and 



(3.19) 



The scattergrams in Fig. and Fig. |T^, based on 40000 generated 3-jet events each, display 
the ideal [i values (ordinate) versus y! (abscissa), for cm. energies of 80 and 200 GeV. 



Obviously the vast majority of events have ^ > . This is also seen in Fig. 16, where the 

A^' = ^ - /i' (3.20) 
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Figure 16: Distribution of A^' = ji — ji' for generated 3-jet events of 40 GeV cm. energy 
of 40 GeV events are plotted. In addition, the scattergrams show a depletion of events 



near ^' = 0.25. In Fig. 16 the mean value is (A^u') = 0.0479, and the standard deviation is 
a = 0.0171. 




Figure 17: A/x' (3.20) (ordinate) vs /i' (3.19) (abscissa) of 3-jet events for different energies. 



In order to obtain a more correct approximation to the value of as a function of the 
measured quantities fi' and the cm. energy E, we parametrized A//' ( 3.2C| ) as a function 
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of these variables. Fig. 17 shows the average values of A/i' for given spans of n' , and for 
different energies. In addition to the span of ^' below the value 0.05, the shaded parts are 
also not used, in order to avoid the //' values from the depleted regions. The straight lines, 
corresponding each to a given energy, were obtained from fitting by "Minuit" |17], yielding 
= 0.19 for 11 degrees of freedom. The results are: 



E + e 



[fi' + a) 



(3.21) 



where the energy E is in GeV, and the numerical values of the parameters are: c = 7. 167 it 
1.204; e = —5.16 it 4.09 and a = 0.1098 ± .0404. As in the case of the parameterization 



(3.16), the extrapolation to energies lower than 40 GeV is not safe. The more correct 
approximation of as a function of the measured values is therefore 



/i2 = ^' + A2 {ii',E) 



(3.22) 



The straight line A^' + fi' = 0.25 from Fig. ^ limits the maximal value of fi2 to the 
physical limit of 0.25, meaning that if one obtains fi2 > 0.25, the value of 0.25 should he 
used. Otherwise, there are no restrictions on the /i' values, for using ( |3.22| ). 
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Figure 18: ^ (ordinate) vs ^2 (abscissa) of Figure 19: ^ (ordinate) vs ^2 (abscissa) of 
3-jet events with 80 GeV cm. energy. 3-jet events with 200 GeV cm. energy. 



After applying the correction of ( |3.22| ) one obtains the scatter plots Fig. |l^ and Fig. 19 



which should be compared with Fig. |lj and Fig. 15, using uncorrected n' values. 
The display of the corrected value 



A/i2 



(3.23) 



presented in Fig. ^ should be compared to the - uncorrected of Fig. 16. After correct- 
ing, one obtains improved values of the mean: —0.00248, and of the standard deviation: 
0.01573. Even so, these values are worse than those obtained using the tensor of rank-2; 
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Figure 20: Distribution of A^2 — fJ, ~ H2 for generated 3-jet events of 40 GeV cm. energy 
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Figure 21: Distribution of ^12/ tJ^i for gener- Figure 22: Mean value and standard devi- 
ated 3-jet events of 80 GeV cm. energy. ation of A'2//ii for different cm. energies. 



for comparison see Fig. ^ and Fig. For this reason the values of ^1 are more reliable than 
those of ^2- Nevertheless, their ratio is close to unity, as expected from ( 3.1^ ), and shown 
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in Fig. 21, for energy of 80 GeV. The mean value in this case is 0.9959 and the standard 
deviation is 0.0592, which makes it very convenient for picking up 3-jet events over a large 
background. This kind of behaviour is characteristic for all the studied energies, from 40 
to 200 GeV, as shown in Fig. [2^. 



3.3 Real data 

This method differs from the generally accepted ones for jet analysis, therefore it is neces- 
sary to discuss how to implement it on real data. The method is applicable for 3-jet events 
only, therefore we start by reviewing the selection criteria for such events. 

The momentum tensors require to have all the particles originating from the event's 
vertex, boosted to their common centre of mass system. An event should be accepted 
only if the centre of mass energy E is above a limit of at least 40GeV, in order to obtain 
accurate enough results. The momentum tensor of rank-2 Qij ( |1.1] ) should be calculated 
and its eigenvalues (A's) obtained and ordered (|1.4|). The value of /xi ( p. 12 ), should be 



calculated and used as an approximation to the idealized fi value (^.l|). By definition its 
value is between and 0.25 . An event should be accepted only if its value is above a limit 
of at least 0.05, in order to avoid 2-jet events. 

The momentum tensor of rank three Qijk ( ^-12 ) should be calculated, and its corre- 



sponding tensor of rank two Rij ( p.l6| - ^T7 ) obtained. The eigenvalues (u) of the latter 



should be ordered, and the value of fx' (|3.19D calculated. This value is infiuenced by the 
hadronization, and needs to be transformed to fi2 ( |3.21[ - 3l2^ ), following the related instruc- 



tions in the text. As a result the distribution of (/ii — fi2) should include a peak centred 
close to the zero, and a background mostly toward the positive values. Also the {^2/ ^^-i) 
distribution should show a peak pointing to a value close to one with a background, mostly 
to the left. These peaks correspond to the 3-jet events, and the second one is expected to 
be sharper, therefore more suitable for the selection. 

The transformation to ^2 just described, corresponds to the hadronization of the event 
generator used in the paper. In order to obtain more realistic results, the parameters "a" 
and "c" from ( p.21| ) could be tuned, in order to move the centre of the peaks to the expected 
values. 

After the 3-jet events are selected, the present study makes possible the direct mea- 
surement the hadronization corrections of the /i value ( |3.11 , p. 13 ), which are of theoretical 



and practical importance. This may improve the choice of a more realistic event generator 
for Monte-Carlo studies. 

Comparisons between the results based on the analysis of 3-jet events using the mo- 
mentum tensors, and between the identified jets by the clustering methods, can be made, 
especially on event by event basis. Of course the identified jets obtained from the cluster- 
ing should be boosted to their centre of mass system first. Their momenta could be used 
for calculating ^ of the idealized jets (|2.6|) , and compared with /ii (|3.12|) from the tensor 
analysis. The relative efficiency of both methods could also be studied. 
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4. Summary and conclusions 

The rank-2 momentum tensor, calculated from the momenta of the jets, in a 3-jet event, 
provides important characteristics of the event. This is not so, if the tensor's calculation 
is based on all of the particles originating from the vertex of the event, without assigning 
them to the jets. The paper shows that for high centre of mass energy (> 40 GeV), these 
characteristics could be retrieved within some calculable error limits. The method of this 
procedure was constructed and tested by the analysis of Monte-Carlo generated events. 

In addition, the momentum tensor of rank-3 is studied, and it is shown that together 
with the tensor of rank-2, one obtains an unique signature of 3-jet events. This part is also 
supported by an analysis of generated events. 

Since the use of the rank-3 momentum tensor, together with this of rank-2, supplies 
signatures of 3-jet events, it is suggestive to study, if tensors of higher ranks do supply 
signatures of events with higher number of jets. 
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A. Symmetric planar jets, tensor of rank 2 



From the definition of an idealized, rotationally symmetric planar event of > 3 jets 
( 2.10[) (Fig. ^ and the definition of the momentum tensor of rank-2 one obtains 



N N 

Qn = jjE cos2 ( 2-1) = _!_ ^ [1 + eos (^)] 

n=l n=l 
N N 

Q22 = i E sin^ (^) = 2F E [1 - cos (^)] 

n=l n=l 

N N 
Q12 = Q21 = ^ E sin (^) cos (2-) = ^ ^ sin {^) 

n=l n=l 



(A.l) 



but the final sums from ([A.l|) can be expressed in a closed form, (see e.g. |16| 1.341) : 



n=l 



(A.2) 



sm(f ) 
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which yields for any > 3 



s{27r)sin(^27r) 



N 

E cos (^) = 
n=l 

= <522 = 5 and Qi2 = Q21 = 



1 







1 = 



1 

4 



(A.3) 



B. Symmetric planar jets, tensor of rank 3 



From the definition of an ideahzed, rotationally symmetric planar event of A^ > 3 jets 
( 2.10p (Fig. and the definition of the momentum tensor of rank-3 (2.12), one obtains 



N 

Qui = F E C0S3 (2^) 
n=l 
N 

Q222 = ^ E sin3 (2Hi) 

n=l 



In the case of A^=3, by use of (|B.lD and (|2.2C1| ), one obtains 



Qm(Ar = 3) = I 
Q222 (A^ = 3) = 



II = 0.25 



as it should for a 3-jet event. 
Since 



cos^ f = ji'i cos if + cos (3(/?)] "I 
sin'^ = 3 [3 sin ip — sin (399)] J 



one obtains for any A^ 



N 



I COS 



' 2-Kn 



) + 



cos 



' 67rn 
^ TV 



)] 



Qui - 4^ E [3 

n=l 

Q222 = 4F E [3sin(2^)-sin(6Hi)] 

n=l 



which yields using ( [A. 21 ), together with the requirement of A^ > 3 : 



Q222 



1 


3 cos TT sin(7r- 










1 


3 sin TT sin(7r- 




~ 4Af 





+ 



" iV , 



i(37r)sin(37r+^) 



1 
4iV 



sin(37r)sin(^) 
4A^sin(f ) 



(B.l) 



(B.2) 



(B.3) 



(B.4) 



(B.5) 



Qui and Q222 vanish for any A^ > 3. For A^ = 3, Qui and Q222 are undefined in (|B7 
and they should be taken from (|B.2|). 
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